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UNCLASSIFIED  ABSTRACT 


THE  BLUNT  BODY  PROBLEM  IN  AN  TU-0059(SGblG-7G)-2 

EXPANDING,  NONUNIFORM  FLOW  FIELD  70  SEP  07 

by  P.  G.  Crowell 


The  method  of  integral  relations  has  been  used  to  compute  the  inviscid  flow 
over  a  blunt  body  immersed  in  a  1.  nuniform,  expanding  flow  Held.  Emphasis  has 
been  placed  on  defining  the  How  environment  experienced  by  a  nose  Up  model  in  the 
exhaust  plume  of  an  underexpanded  arc  jet  or  rochet  engine.  Accordingly,  the  integral 
method  has  been  formulated  for  an  arbitrary  description  of  a  free  stream  which 
possesses  both  axial  and  radial  gradients  in  .Mach  No.,  flow  angle,  total  pressure,  and 
total  enthalpy.  , 

Computations  are  presented  for  two  different  descriptions  of  the  free  stream 
(spherical  source  flow  and  the  exhaust  plume  from  a  contoured,  Mach  2  nozzle)  and 
two  body  shapes  (sphere  and  circular  disk).  Results  indicate  a  first  order  effect  on 
all  the  blunt  body  characteristics. 

The  influence  of  the  free  stream  derivatives  of  Mach  No.  and  flow  angle  has  been 
examined  and  found  to  exert  only  a  second  order  effect  upon  the  solutions,  thereby 
allowing  a  considerably  simplified  description  of  the  free  stream.  Move  over,  it  lias 
been  demonstrated  that  the  influence  of  an  expanding  flow  lie  1  cl  upon  a  blunt  body  may 
be  approximately  characterized  in  terms  of  the  flow  ancle  gradient  normal  to  the 
plume  axis.  Consequently,  definition  of  the  plume  flow  field  is  required  only  along 
its  axis. 

The  possibility  of  sonic  point  movement  off  a  sharp  corner  onto  the  face  of  a 
circular. disk  has  been  investigated  and  found  to  be  entirely  possible  for  a  nonuniform 
free  stream.  To  identify  the  magnitude  of  the  free  stream  flow  divergence  required 
to  move  the  sonic  point  off  the  corner,  a  set  of  solutions  was  obtained  for  an  infinite 
circular  disk  in  a  spherical  source  flow. 

Model  size  constraints,  for  the  simulation  of  uniform  flow  pressure  and  heating 
distributions,  have  been  identified  for  ablation  testing  in  exhaust  plumes.  Furthermore, 
the  region  of  a  plume  flow  field  immediately  downstream  of  the  Mach  rhombus  has  been 
identified  as  the  worst  possible  test  location  if  simulation  of  a  uniform  flow  environ¬ 
ment  is  desired. 

Requirements  have  been  identified  for  additional  cold  flow  experiments  to 
define  the  capabilities  and  limitations  of  simulating  uniform  flow  environments  in  under- 
expanded  exhaust  plumes.  (Unclassified  Report) 
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SECTION  1 


INTRODUCTION 


There  are  several  fluid  mechanics  problems  of  current  interest  where  it  is  necessary 
to  compute  the  flow  about  a  blunt  nosed  body  in  a  nonuniform  free  stream.  These 
problems  may  be  considered  to  fall  into  two  different  categories,  those  where  the  free 
stream  is  essentially  parallel  and  where  the  free  stream  is  diverging.  Nonuniform, 
parallel  flow  is  encountered  by  a  body  when  it  is  immersed  in  a  constant  pressure 
shear  layer  such  as  a  wake  or  jet  flow.  A  typical  example  of  interest  consists  of 
deploying  objects  into  the  wake  of  a  reentry  vehicle. 

When  a  body  is  immersed  in  an  expanding  flow  field,  such  as  a  conical  nozzle  or 
exhaust  plume,  the  free  stream  contains  gradients  in  both  the  flow  angle  and  gas 
properties.  Several  investigations  have  been  made  concerning  the  effects  of  a  small 
amount  of  flow  divergence  upon  the  flow  field  properties  of  a  body  in  a  wind  tunnel. 
Reference  1,  for  example,  has  used  an  inverse  method  with  the  method  of  character¬ 
istics,  to  compute  the  inviscid  flow  over  spherical  and  elliptical  bodies  in  a  spherical 
source  flow.  Attention  was  restricted,  however,  to  cases  where  the  flow  divergence 
was  only  a  few  degrees. 

The  interest  in  high  performance  reentry  vehicles  has  led  to  considerable  testing  of 
ablative  and  transpiration  nosetips  in  highly  underexpanded  free  jets.  A  rapidly 
expanding  free  jet  possesses  large  gradients  in  flow  properties  and  flow  angle,  except 
in  the  Mach  rhombus  formed. at  the  exit  of  a  contoured  nozzle.  It  is  appai'ent,  there¬ 
fore,  that  in  many  ground  test  facilities,  the  flow  over  the  nosetip  will  possess  large 
gradients  in  flow  properties  in  both  the  radial  and  axial  directions.  Depending  upon 
the  relative  sizes  of  the  nozzle  exit  diameter  and  model  diameter  and  the  distance  of 
the  model  from  the  exit  plane,  free  stream  flow  divergence  angles  of  up  to  twenty 
degrees  may  be  encountered. 
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In  spite  of  the  interest  in  such  testing,  there  does  not  appear  to  be  any  systematic 
investigation  of  the  effects  of  large  free  stream  flow  divergence  angles  upon  the  blunt 
body  problem.  Ref  (2)  used  the  method  of  integral  relations  coupled  to  a  character¬ 
istics  solution  to  compute  the  flow  field  for  a  sphere  in  a  rocket  exhaust,  and  Ref  CJ) 
has  used  a  time  dependent  technique  for  a  concave  body  in  an  exhaust  plume.  However, 
in  each  case,  only  one  or  two  calculations  were  carried  out  for  a  specific  body  con¬ 
figuration  and  test  geometry. 

It  is  the  purpose  of  this  report  to  investigate  the  flow  of  a  nonuniform,  highly  diverging 
free  stream  over  a  blunt  nosed  body.  In  particular,  model  size  constraints  will  be 
identified  for  ablation  testing  in  a  free  jet  and  the  effect  of  nonuniform  flow  upon  the 
simulated  flight  environment  will  be  examined. 


e 
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SECTION  2 


DEFINITION  OF  FREE  STREAM  FLOW  FIELD 


As  mentioned  previously,  the  flow  field  is  to  be  calculated  over  a  blunt  nosed  body 
immersed  in  the  exhaust  plume  of  an  underexpanded  nozzle.  Since  the  nozzle  flow  is 
supersonic,  the  method  of  characteristics  has  been  used  to  compute  the  exhaust 
plume  flow  field  for  a  contoured  nozzle  designed  to  produce  a  parallel  Mach  2  flov. 
at  the  exit  (Ref  13).  Mach  No.  and  flow  angle  distributions  and  their  gradients  in  the 
radial  and  axial  directions  are  plotted  in  Figures  2  through  9.  Figure  1  defines  the 
geometry  for  the  blunt  body  and  the  external  flow  fiel  i. 

In  addition  to  the  free  jet  flow  fields,  computations  have  also  been  made  for  a  spherical 
source  flow.  Mach  No.  and  flow  angle  distributions  for  the  spherical  source  flow  may 
be  described  by  the  following  expressions: 


M 


r  +  8Cos0 

Tan  4*  "  [  -  (5  Sin  0  -  Z) 
o 


At  any  point  on  the  shock  wave; 


It  is  noted  that  the  expression  for  the  Mach  No.  variation  assumes  that  the  velocity  is 
constant  throughout  the  flow  field,  such  that: 


In  the  above  expressions,  the  subscript  o  refers  to  source  flow  properties  which 
would  exist  at  the  location  of  the  body  stagnation  point  if  the  body  were  not  immersed 
in  the  flow  field. 
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The  calculation  of  the  external  flow  field  is  uncoupled  from  the  calculation  over  a  body 
immersed  in  that  flow  field.  Therefore,  it  is  required  that  the  body  shock  wave  not 
influence  the  free  stream  in  a  region  where  the  free  stream  definition  is  required  for 
continuing  the  computation  of  the  shock.  This  may  be  thought  of  in  terms  of  a  restric¬ 
tion  on  the  free  stream  velocity  and  flow  divergence  angle,  such  that  the  component  of 
the  free  stream  velocity  vector,  normal  to  the  shock  wave,  always  remain  supersonic. 
That  is,  it  is  required  that: 

Mffl  Sin (/3-<f>)  >  1 

It  is  possible,  therefore,  that  there  may  exist  some  combination  of  body  shapes  and 
free  stream  flow  conditions  which  will  preclude  uncoupling  the  body  and  free  stream 
flow  field  calculations  beyond  some  limiting  distance  downstream  of  the  sonic  point. 

In  fact,  there  exists  some  limiting  flow  divergence  angle  such  that  the  above  restric¬ 
tion  wall  be  violated  before  the  subsonic  region  of  the  flow  field  is  closed.  In  this 
event,  the  two  flow  fields  wall  be  completely  coupled  and  the  computation  will  fail. 
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SECTION  3 


METHOD  OF  INTEGRAL  RELATIONS  FOR  A  NONUNIFORM  FREE  STREAM 


The  method  of  integral  relations  (one  strip  approximation)  has  been  used  to  make  the 
computations  in  this  report.  It  was  decided  not  to  use  the  inverse  method  of  Ref  1 
because  of  the  difficulty  of  prescribing  the  shock  shape  for  a  nonuniform  free  stream, 
particularly  where  the  external  flow  divergence  angle  is  large  enough  to  cause  inflec¬ 
tions  in  the  shock  shape.  Furthermore,  Inouyc  (Ref  1)  has  reported  encountering 
numerical  difficulties  when  the  flow  divergence  increased  beyond  a  few  degrees.  In 
addition,  the  inverse  method  cannot  be  used  to  handle  a  body  with  a  sharp  sonic  corner. 
For  these  reasons,  the  direct  method  has  been  used. 

In  the  development  of  the  method  of  integral  relations  which  follows,  the  fluid  has  been 
assumed  to  be  a  perfect  gas  where  the  ratio  of  specific  heats  is  constant  throughout  the 
flow  field.  The  free  stream  possesses  both  axial  and  radial  gradients  in  Mach  No. ,  flow 
angle,  total  enthalpy,  and  total  pressure.  In  the  following  discussion,  a  general  outline 
is  presented  that  emphasizes  the  differences  between  uniform  and  nonuniform  flow. 

A  detailed  description  is  included  in  Appendix  1. 

The  continuity  and  normal  momentum  equations  are  written  in  body  coordinates  and 
integrated  along  a  direction  normal  to  the  body  surface  out  to  the  shock  wave. 

Assuming  that  the  resultant  integrands  are  linear  in  the  normal  direction,  n  ,  and 
using  the  notation  of  South  (Ref  4),  we  obtain: 
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From  geometrical  considerations  (see  Figure  1),  the  third  equation  is  obtained: 
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It  is  noted  that  equations  1,  2,  and  3  are  valid  for  any  description  of  the  free  stream. 

It  is  apparent  that  the  free  stream  properties  enter  the  differential  equations  through 
the  use  of  the  oblique  shock  relations  to  define  the  quantities  behind  the  shock.  There¬ 
fore,  to  evaluate  the  quantities  behind  the  shock  (subscript  1)  the  shock  angle  in  the 
oblique  shock  equations  is  simply  /3  -  (see  Figure  1).  The  only  substantial  eficct 
of  the  nonuniform  free  stream  upon  the  formulation  is  the  presence  of  additional  terms 
in  the  differential  equations  (Ref  6).  These  additional  terms  enter  the  equations  when 
the  derivatives  of  p^u^  Vj  and  are  evaluated.  For  example,  consider  P1U2V1: 
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From  the  oblique  shock  relations,  it  is  seen  that:  Fj  ~  F^  ^/3,  6 ,  <f>, 


G 
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It  is  apparent,  therefore,  that  differentiation  of  F1  and  will  yield  derivatives 

of  the  flow  angle,  Mach  No. ,  and  total  pressure  which  would  be  zero  for  a  uniform  flow. 
Evaluation  of  ^  (t^u  )  w*^  yield  similar  terms.  The  complete  details  are  contain'd 
in  Appendix  1. 

The  presence  of  these  additional  terms  in  the  differential  equations  complicates  the 
description  of  the  free  stream  environment.  If  the  external  flow  is  described  by 
analytic  expressions  for  Mach  No.,  flow  angle,  total  pressure,  and  enthalpy,  then  these 
expressions  may  be  differentiated  to  obtain  the  additional  gradients.  If,  however,  the 
free  stream  definition  is  obtained  by  either  curve  fitting  or  table  look  up  and  interpola¬ 
tion  of  a  method  of  characteristics  solution,  the  requirement  for  evaluating  these  addi¬ 
tional  terms  involves  a  considerable  effort.  It  will  be  of  interest,  therefore,  to  assess 
the  effect  of  ignoring  the  derivatives  of  the  free  stream  properties.  This  will  be 
discussed  further  in  Section  4. 

As  in  the  uniform  flow  case,  Eq  1  and  2  yield  indeterminate  expressions  for  the  initial 
derivatives  of  j3  and  uQ  .  Application  of  L'hopitals  ride  yields  the  following  expres¬ 
sions  for  a  nonuniform  external  flow:  (See  Appendix  2. ) 


Is] 

1  dP 

-  g  P0"P1  _ 

P, 0 

pl> 

ds 

ds  p  v2 

co  va? 

1 

duo 

Pj(l  +  K8)  (  P0- 

P1 

ds 

V  [  p. 

V2 

V  CO 

It  is  quite  intei'esting  to  note  that  the  stagnation  point  velocity  gradient  has  exactly  the 
same  form  for  either  uniform  or  nonimiform  flow.  This  is  of  considerable  importance, 
since  it  implies  that  the  stagnation  point  velocity  gradient  may  still  be  obtained  from  a 
measurement  of  the  shock  standoff  distance,  even  in  a  highly  diverging  flow  field. 
Furthermore,  although  the  method  of  integral  relations  does  not  accurately  predict  the 


stagnation  point  velocity  gradient  (Ref  7),  it  may  be  obtained  from  the  uniform  flow 
correlation  of  Kaattari  (Ref  8),  even  for  a  nonuniform  flow.  Kaattari's  correlation 
requires  the  shock  stand  off  distance,  but  this  may  be  obtained  accurately  from  the 
method  of  integral  relations. 


The  initial  derivative  for  the  shock  inclination  angle,  /? ,  differs  from  the  uniform 
flow  result  by  the  additional  term  for  the  derivative  of  the  flow  divergence  angle,  <f> . 
Depending  upon  the  magnitude  of  4>  and  its  gradient,  the  external  flow  divergence  may 
have  a  significant  effect  upon  boundary  layer  edge  properties  and,  thus,  the  heat  trans¬ 
fer.  This  effect  is  due  to  the  fact  that  entropy'  swallowing  is  dependent  upon  both  the 
shock  angle  and  its  gradient  normal  to  the  flow  field  axis.  It  is  possible,  of  course, 
that  this  effect  may  be  overwhelmed  by  the  presence  of  a  large  entropy  gradient  in  the 
free  stream  (i.c.,  nonconstant  total  enthalpy  and  total  pressure). 

In  order  to  evaluate  the  derivatives  of  the  free  stream  properties  with  respect  to  the 
s  direction,  a  transformation  is  required  between  the  body  centered  coordinate  sys¬ 
tem  and  the  coordinate  system  of  the  external  flow.  The  free  stream  properties  shall 
be  specified  in  either  a  cylindrical  or  spherical  coordinate  system.  In  all  cases, 
circumferential  symmetry  is  assumed.  The  gradient  of  any  scalar  property  P  may¬ 
be  written  in  cylindrical,  spherical,  or  body  centered  coordinates  as  follows: 


Cylindrical 

VP  = 

r* 

dr  r 

dp  — 
a  z  ez 

Spherical 

VP  = 

dP  — 
d( 

1  dP  — 

f  a*  c* 

Body 

VP  = 

1 

dp  Trt  A  dP  — ► 

1  +  Kn 

ds  s  dn  en 

The  dot  products  of  the  various  unit  vectors  are  as  follows: 
•  c^  =  Sin  8  e"^  •  =  Cos  9 

eg  •  =  Cos  (9-<f>)  c^  •  =  Sin  (9-<f>) 


t 
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Taking  the  appropriate  dot  products  between  the  gradients  and  the  unit  vectors,  the 
transformations  are  obtained. 


For  a  spherical  external  flow: 


-  (l'1  Kn)  Cos  ( 8-<t> )  — j-  + 


Sin  -|| 


For  a  cylindrical  external  flow: 

=  (1  +  Kn)  Sin  0  ~  +  (1  +  Kn)  Cos  d  —  (6) 

For  a  spherical  source  flow,  the  gradients  in  Mach  No,  and  flow  angle  take  on  particu¬ 
larly  simple  forms  because  the  thermodynamic  state  and  flux  properties  are  independent 
of  <t>  and  the  flow  angle  is  independent  of  f  . 

dM  dM 

=  (1  +Kn)  Cos 


=  (1  +  Kn) 


Sin  (d-4>) 


Eq  1-3,  along  with  the  boundary  conditions  (Eq  10)  constitute  a  two  point  boundary  value 
problem  whose  nature  is  unchanged  by  the  fact  that  the  external  flow  field  is  non- 
uniform.  Therefore,  as  in  the  uniform  flow  case,  the  solution  is  obtained  by  iterating 
upon  the  initial  value  of  the  shock  standoff  distance  until  the  correct  sonic  point 
behavior  is  achieved.  A  body  whose  sonic  point  is  located  at  a  sharp  corner  exhibits  a 
singularity  in  the  velocity  gradient  at  that  point  and  the  correct  behavior  for  the  solu¬ 
tion  is  simply  that  the  sonic  velocity  be  attained  at  the  corner.  For  a  body  with  sonic 
point  located  on  a  smooth  contour,  it  is  required  that  the  velocity  gradient  remain 
finite  as  the  sonic  velocity  is  approached.  For  a  comprehensive  discussion  of  the 
behavior  of  the  solution  in  the  vicinity  of  the  sonic  point,  Ref  10  may  be  consulted. 


uQ(0)  =  0 


V3  *>  =  u0 


/3(0)  =  it/ 2 
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In  light  of  the  fact  that  the  computational  aspects  of  the  problem  are  unchanged  in  the 
presence  of  a  nonuniform  free  stream,  it  was  decided  to  modify  existing  computer 
programs  to  account  for  nonuniformities  in  the  external  flow.  The  programs  that  were 
used  in  the  study  were  written  by  Jerry  C.  South  of  NASA,  Langley.  The  program  to 
compute  the  flow  field  for  bodies  with  sonic  corners,  in  a  uniform  free  stream,  is 
documented  in  Ref  4.  The  second  computer  program,  used  to  obtain  the  sonic  saddle 
point  convergence  for  smooth  bodies,  has  been  informally  documented  by  South  (Ref  o» 
Both  programs  used  a  fourth  order  Runge  Kutta  integration  scheme  with  a  variable 
step  size. 
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SECTION  4 


SPHERICAL  SOURCE  FLOW  RESULTS 


It  is  well  known  that  the  plume  from  an  underexpancled  nozzle  may  be  simulated  with 
reasonable  accuracy  by  a  spherical  source  flow  located  at  the  virtual  origin  of  the 
plume  flow  field.  Furthermore,  in  the  case  of  a  spherical  source  flow,  the  effect  of 
the  diverging  free  stream  impinging  upon  a  blunt  body  may  be  characterized  by  a  .single- 
parameter,  the  ratio  of  the  characteristic  body  length  scale  to  the  source  distance 
from  the  body. 

For  these  reasons,  an  extensive  set  of  calculations  has  been  carried  out  for  a  sphere 
and  a  circular  disk  immersed  in  a  spherical  source  flow.  The  calculations  covered  a 
Mach  No.  range  from  2.0  to  20.  0  and  %ralues  of  the  specific  heat  ratio  of  1. 2,  1.4  and 
1.6667. 

Before  presenting  the  results  of  the  calculations,  a  qualitative  description  of  the  effects 
of  a  diverging  free  stream  upon  the  blunt  body  flow  field  is  in  order.  The  critical 
difference  between  a  blunt  body  flow  field  with  a  parallel  external  flow  and  one  with  a 
diverging  flow  is  that  the  shock  angle,  (i.e.,  0  -  <f>)  decreases  faster  for  the  diverging 
flow,  for  the  obvious  reason  that  4>  is  non  zero.  Therefore,  the  pressure  behind  the 
shock,  and  thus  also  on  the  body  surface,  decreases  faster  than  for  the  parallel  flow 
case.  This  leads  to  an  increase  in  the  surface  velocity  gradient  which  causes  a  move¬ 
ment  of  the  sonic  point  towards  the  stagnation  point  (except  for  a  body  whose  sonic 
point  is  fixed  at  a  sharp  corner). 

The  behavior  of  the  shock  standoff  distance  may  be  inferred  from  one-dimensional 
streamtube  continuity  considerations.  Since  the  sonic  point  location  on  the  shock 
moves  closer  to  the  axis  as  the  flow  divergence  increases,  the  free  stream  flow  angle 
always  remains  relatively  small  up  to  the  sonic  point.  Therefore,  the  mass  flux  pass¬ 
ing  through  the  shock  at  any  radial  distance  from  the  axis  is  approximately  the  same 
for  either  the  uniform  or  non-uniform  flow,  provided  the  free  stream  Mach  numbers  are 
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not  significantly  different.  Since  the  mass  flux  per  unit  area  between  the  shock  wave 
and  body  surface  is  greater  for  the  nonuniform  flow,  continuity  of  mass  requires  that 
the  flow  area  (and  thus  the  shock  standoff  distance)  be  smaller  than  for  the  correspon¬ 
ding  uniform  flow. 

To  summarize,  the  effect  of  a  diverging  free  stream  is  to  increase  the  stagnation 
point  velocity  gradient,  decrease  the  surface  pressure,  move  the  sonic  point  closer  to 
the  axis,  and  decrease  the  shock  standoff  distance. 

Typical  results  for  the  shock  shape,  pressure  and  heat  flux  distributions  on  a  spherical 

body  are  presented  in  Figures  10  through  13.  The  laminar  heat  flux  is  Lee's  cold  wall 

distribution  and  the  turbulent  heating  is  Vaglio  Laurin's  formulation.  It  is  apparent 

that  for  values  of  the  source  flow  parameter  y-  of  up  to  0.5,  the  diverging  free  stream 

*  o 

has  a  first  order  effect  on  all  the  blunt  body  properties.  In  Figures  14  and  15,  the 
sonic  point  location  and  shock  standoff  distance  are  plotted  for  a  range  of  Mach  No. 

The  shock  standoff  distance  and  stagnation  point  velocity  gradient  may  vary  by  almost 
a  factor  of  two  over  their  uniform  flow  values.  By  taking  the  square  root  of  the 
velocity  gradient  ratio  in  Figures  IS  and  19,  an  estimate  of  the  increase  in  the  laminar 
stagnation  point  heat  flax  may  be  obtained.  This  indicates  that  for  low  supersonic 
Mach  values,  the  stagnation  point  heating  may  be  increased  by  thirty  five  percent  over 
the  corresponding  uniform  flow  value.  The  sensitivity  to  Mach  No.  at  low  supersonic 
values,  which  is  indicated  in  Figures  1G  and  18  is  due  to  the  sensitivity  of  the  normal 
shock  density  ratio  to  low  supersonic  values  of  the  incident  Mach  No.  Table  1  contains 
a  summary  of  the  solutions  for  a  sphere  in  a  spherical  source  flow. 

Since  the  results  of  Ref  1  were  obtained  bj'  an  exact  numerical  solution  of  the  inviscid 
equations,  it  is  of  interest  to  compare  them  with  the  results  of  the  integral  method  of 
this  report.  The  comparison  of  shock  standoff  distance,  stagnation  point  velocity 
gradient,  and  pressure  distribution  is  contained  in  Figures  20  through  22.  The 
standoff  distance  and  pressure  distributions  are  in  excellent  agreement.  The  stagna¬ 
tion  point  velocity  gradient  is  in  error  by  about  18  percent,  as  expected  (Ref  7). 
However,  the  values  of  the  stagnation  point  velocity  gradient  normalized  with  the 
uniform  flow  value  are  in  excellent  agreement.  It  is  felt,  therefore,  that  the  velocity 
gradient  ratio  may  lie  accurately  calculated  by  the  integral  method. 
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The  use  of  Newtonian  impact  theory  to  compute  pressure  distributions  on  blunt  bodies 
is  a  simple,  yet  accurate  technique  for  uniform  flows.  Newtonian  theory  might  also  be 
expected  to  yield  adequate  predictions  for  bodies  in  an  expanding  flow  field.  Using  the 
source  density  and  flow  angle  along  the  body  surface,  the  following  expression  is 
obtained  for  the  pressure  distribution  on  a  sphere  in  a  spherical  source  flow. 


STAG 


Cos“0 


i-4>)  o 

\  R 


,  Tan  4>  =  —  Cos  6 


The  comparison  of  the  pressure  distributions  is  contained  in  Figure  23.  It  is  seen  that 
the  Newtonian  prediction  is  in  reasonable  agreement  (within  10  percent)  with  the 
results  of  the  integral  method,  except  near  the  shoulder  of  the  sphere.  However,  the 
integral  method  itself  is  inaccurate  in  this  region.  This  is  clearly  evident  from  the 
uniform  flow  calculation  where  the  integral  method  yields  pressures  which  are  some¬ 
what  higher  than  would  be  obtained  from  a  characteristics  calculation.  Therefore,  for 
small  values  of  6 ,  the  Newtonian  result  may  be  more  accurate  than  the  integral 
method. 

As  was  mentioned  earlier,  one  of  the  disadvantages  of  the  integral  method  is  that  the 
free  stream  derivatives  of  Mach  No.  and  flow  angle  must  be  specified  as  well  as  the 
Mach  No.  and  flow  angle  themselves.  Since  the  most  significant  aspect  of  the  diverging 
free  stream  is  its  effect  on  the  shock  angle,  it  might  be  expected  thaf  these  derivatives 
would  exert  only  a  second  order  effect  on  the  solution.  In  order  to  ascertain  the 
importance  of  these  terms  upon  the  solution,  a  set  of  solutions  was  obtained  where  the 
derivatives  of  free  stream  Mach  No.  and  flow  angle  were  ignored.  The  results  for  the 
shock  standoff  distance  and  stagnation  point  velocity'  gradient  are  presented  in  Fig¬ 
ures  24  through  25.  The  values  are  normalized  with  the  corresponding  exact  values. 
Results  indicate  that  ignoring  the  free  stream  gradient  terms  underpredicts  the  stag¬ 
nation  point  velocity  gradient  and  overpredicts  the  shock  standoff  distance.  The  error 
increases  with  increasing  values  of  the  source  flow’  parameter  and  is  approximately 

ten  to  fifteen  percent  for  —  =  0.  5.  The  effect  on  pressure  distribution  is  of  the  same 
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order  and  is  indicated  by  the  dashed  line  in  Figure  10.  It  is  felt,  therefore,  that  in 
situations  when-  it  is  difficult  to  evaluate  the  free  stream  derivatives  that  they  may  be 
ignored  without  introducing  large  errors  into  the  results. 

In  light  of  the  substantial  change  in  the  sonic  point  location  for  a  sphere  in  a  diverging 
flow  field  (see  Figure  11),  it  becomes  pertinent  to  examine  under  what  conditions  the 
sonic  point  will  move  off  a  sharp  corner  onto  the  face  of  a  blunt  body.  For  a  body  sur¬ 
face  with  non  zero  curvature,  the  sonic  point  location  is  computed  as  previously  dis¬ 
cussed  and  truncating  the  body  forward  or  aft  of  this  point  determines  whether  the 
sonic  point  is  located  at  a  corner  or  on  a  smooth  contour. 

The  situation  for  a  body  with  zero  surface  curvature,  (i.e.,  a  sharp  cone  or  circular 
disk)  is  not  quite  so  obvious.  The  contention  that  the  sonic  point  must  always  remain 
at  a  sharp  corner,  for  a  body  with  zero  surface  curvature,  is  based  on  the  argument 
that  if  the  sonic  point  motes  forward  of  the  corner  then  the  corner  can  no  longer 
influence  the  subsonic  flow  field.  In  other  words,  there  no  longer  exists  any  length 
scale  to  the  problem.  That  this  is  physically  reasonable  is  obvious  for  a  uniform 
free  stream. 

However,  for  a  nonuniform  free  stream,  this  argument  is  no  longer  valid  since  there 
are  length  scales  associated  with  the  gradients  in  the  external  flow.  For  a  spherical 
source  flow’  or  an  exhaust  plume,  the  obvious  length  scale  is  the  distance  of  the  body 
from  the  source  or  nozzle  exit  plane.  For  a  flow  which  is  nonuniform  but  parallel, 
there  exists  some  length  scale  associated  with  the  Mach  No.  gradient. 

In  order  to  examine  the  behavior  of  a  blunt  body  with  a  sonic  corner,  in  an  expanding 
flow  field,  a  set  of  solutions  has  been  generated  for  a  circular  disk  in  a  spherical 
source  flow.  A  typical  set  of  results  for  shock  shape,  pressure,  and  heat  flux  distri¬ 
butions  is  contained  in  Figures  26  through  28.  It  is  noted  that  the  shock  standoff  dis¬ 
tance  at  the  sonic  point  is  only  weakly  dependent  on  the  nonuniform  free  stream.  This 
is  due  to  the  fact  that  the  mass  flux  and  the  mass  flux  per  unit  area  through  the  "throat" 
arc  approximately  the  same  for  both  uniform  and  nonuniform  flow.  This  is  a  conse¬ 
quence  of  body  geometry  rather  than  external  flow  dictating  the  location  of  the  sonic 
point.  Therefore,  the  "throat”  area  and  the  standoff  distance  are  virtually  unaffected 
by  the  diverging  free  stream. 
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The  results  for  the  shock  standoff  distance  and  stagnation  point  velocity  gradient  and 
their  sensitivity  to  Mach  No.  and  specific  heat  ratio  are  contained  in  Figures  29  through 
33.  Table  2  contains  a  summary  of  the  solutions  for  a  circular  disk  in  a  spherical 
source  flow. 

It  should  be  mentioned  that  the  scaling  between  the  shock  standoff  distance  and  the  dis¬ 
tance  ‘o  the  sonic  point,  which  exists  for  bodies  with  zero  curvature  in  a  uniform  flow, 
is  no  longer  valid  for  nonuniform  external  flow.  In  Ref  10,  it  is  pointed  out  that  the 
shock  standoff  distance  for  a  body  with  zero  surface  curvature  is  linearly  proportional 
to  the  distance  to  the  sonic  point.  Therefore,  only  one  integration  out  to  the  sonic 
point  is  required  to  establish  the  proportionality  constant.  However,  since  there  exists 
an  additional  length  scale1  associated  with  the  nonuniform  free  stream,  the  solutions 
must  be  obtained  by  iterating  upon  the  shock  standoff  distance,  just  as  for  bodies  with 
finite  curvature. 

It  was  mentioned  earlier  that  Kaattari's  correlation  (Ref  8)  for  the  stagnation  point 
velocity  gradient  should  be  valid  for  nonuniform  flows,  since,  at  the  axis,  the  differen¬ 
tial  equation  for  the  velocity  gradient  has  the  same  form  irrespective  of  the  nature  of 
the  external  flow.  This  is  seen  to  be  the  ease  in  Figure  34,  where  the  velocity 
gradient  ratios  for  a  sphere  and  circular  disk  are  compared  with  the  correlation. 

As  discussed  above,  it  is  desired  to  determine  the  value  of  the  source  flow  parameter 
for  which  the  sonic  point  will  move  off  the  corner  of  the  circular  disk  and  onto  its 
face.  When  such  movement  occurs,  as  far  as  the  subsonic  portion  of  the  flow  field  is 
concerned,  the  disk  is  infinite  in  extent.  The  problem,  therefore,  is  that  of  an  infinite 
circular  disk  in  a  spherical  source  flow  with  the  source  distance  as  the  characteristic 
length.  The  sonic  point  convergence  for  this  case  is  exactly  the  same  as  for  any 
smooth  body,  namely,  that  the  sonic  point  velocity  gradient  remain  finite  as  the  sonic 
velocity  is  approached. 

The  results  for  the  infinite  circular  disk  are  contained  in  Figures  35  through  39.  It  is 
noted  that  since  the  disk  is  infinite  in  extent,  the  shock  curvature  is  away  from  the 
body.  For  a  finite  disk  whose  sonic  point  has  moved  onto  the  face  of  the  disk,  the 
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shock  will  bo  concave  in  the  vicinity  of  the  stagnation  point  and  convex  as  the  corner 
is  approached.  Actually,  an  examination  of  the  shock  shapes  in  Figure  28  indicates 
that  the  shock  is  becoming  concave  in  the  vicinity  of  the  stagnation  point  for  valu*  s  of 
the  source  flow  parameter  of  0.333.  The  sonic  point  location  on  the  infinite  circular 
disk  is  shown  in  Figure  35.  Therefore,  for  values  of  the  source  flow  parameter 
greater  than  the  sonic  distance  for  the  infinite  circular  disk,  the  sonic  point  will  move 
off  the  sharp  corner.  Table  3  contains  the  stagnation  and  sonic  point  properties  for  the 
solutions. 

Although  the  inviscid  flow  equations  admit  solutions  with  the  sonic  point  located  ahead 
of  the  sharp  corner  on  a  circular  disk,  there  is  some  question  as  to  whether  viscous 
effects  will  render  the  solution  unstable.  The  available  experimental  data  (Ref  11  and 
12)  for  an  undcrexpanded  jet  impinging  upon  a  flat  plate  appear  to  be  limited  to  sur¬ 
faces  which  are  infinite  in  extent,  with  respect  to  the  exhaust  plume.  Additional 
experiments  with  a  finite  plate  are  required  to  answ'er  the  stability  question. 
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SECTION  5 


EXHAUST  PLUME  CALCULATIONS 


In  order  to  obtain  some  results  which  would  be  of  quantitative  use  in  defining  the 
environment  on  a  nose  tip  in  an  arc  jet  flow  field,  solutions  have  been  generated  for  a 
sphere  in  the  exhaust  plume  of  a  contoured  nozzle.  The  nozzle  was  assumed  to  deliver 
a  uniform,  parallel,  Mach  2  flow  at  its  exit  plane.  Total  enthalpy  and  pressure  were 
assumed  to  be  constant.  The  exhaust  was  assumed  to  be  a  perfect  gas  (  7  --  1.10)  with 
a  ratio  of  chamber  pressure  to  ambient  pressure  of  100. 

The  method  of  characteristics  (Ref  13)  was  used  to  compute  the  exhaust  plume  flow 
field  and  the  results  are  contained  in  Figures  2  through  9.  In  Figure  2,  the  physical 
size  of  the  plume  and  the  location  of  the  boundary  shock  are  shown.  Since  the  flow  is 
uniform  and  parallel  in  the  nozzle  exit  plane,  it  wall  also  be  uniform  inside  a  region 
bounded  by  the  axis  and  the  lead  ray  in  the  Prandtl  Meyer  fan  emanating  from  the 
nozzle  lip.  This  uniform  flow  region  is  indicated  by  the  dashed  line  in  Figure  2. 

The  characteristics  solution  was  written  on  magnetic  tape  and  used  with  a  sur'ace 
interpolation  procedure  to  provide  the  definition  of  the  free  stream  required  for  me 
blunt  body  solution.  The  surface  interpolation  scheme  yielded  a  random  error  of  one 
to  two  percent  for  all  the  How  properties.  This  error,  in  turn,  cause  d  a  random 
error  in  the  surface  velocity  gradient  on  the  order  of  ten  percent.  In  order  to  obtain 
well-behaved  solutions,  the  following  procedure  was  used.  First,  the  shock  wave  loca¬ 
tion  was  approximately  located  by  using  the  interpolated  results  for  the  free  stream 
properties.  Then  the  interpolated  results  were  plotted  and  smoothed  manually.  The 
free  stream  variables  were  then  fit  with  polynomial  expressions  and  the  solution 
was  rerun. 

Fortunately,  the  shock  wave  location  was  not  sensitive  to  the  error  in  the  free  stream 
properties,  so  this  tedious  procedure  only  had  to  be  repealed  once  for  each  case.  In 
any  situation  where  it  was  desired  to  generate  a  large  number  of  solutions,  some  auto¬ 
mated  procedure  would  have  to  be  developed  to  smooth  the  interpolated  values  during 
the  solution. 


( 
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The  computations  for  the  plume  fhw  tu  M  m  in  terms  of  two  parameters,  the  distance 

of  the  sphere  from  the  nuzzle  exit  pi  nu  (~  )  and  the  size  of  the  sphere  relative  to 

/  H  \  '  1  e 

the  nozzle  exit  radius  (  — -  )  .  The  corresponding  uniform  flow  case  occurs  as  the 
sphere  radius  shrinks  to  zero.  A  complete  set  of  results  is  presented  in  Figures  10 
through  12  for  a  sphere  located  one  nozzle  exit  radius  downstream  of  the  nozzle.  At 
this  location,  part  of  the  sphere  lies  inside  the  Mach  rhombus.  Therefore,  there  is  a 
discontinuity  in  the  flow  environment  as  the  bow  shock  cuts  across  the  Prandtl  Moyer 
ray  which  bounds  the  Mach  rhombus.  This  is  reflected  in  the  pressure  distribution 
in  Figure  40.  It  is  apparent  from  the  heat  transfer  and  pressure  that  it  is  not  possible 
to  obtain  uniform  flow  distributions  for  a  sphere  which  is  any  larger  than  half  the 
nozzle  exit  radius.  If  it  is  desired  to  test  larger  bodies  and  still  obtain  heat  transfer 
and  pressure  distributions  which  are  representative  of  those  encountered  for  a  uniform 
fi’ec  stream,  the  body  must  be  placed  further  downstream  of  the  nozzle  exit  plane. 
Accordingly,  a  set  of  pressure  distributions  is  presented  in  Figure  43  for  a  sphere 
which  is  four  nozzle  radii  downstream.  Results  indicate  that  the  sensitivity  to  body 
size  is  substantially  less  at  larger  distances  from  the  nozzle.  This  is  due  to  the 
decrease  in  flow  angle  with  increasing  distance  from  the  nozzle,  for  the  same  radial 
distance  from  the  axis. 

Unfortunately,  since  the  centerline  Mach  No.  increases  quite  rapidly  downstream  of 
the  nozzle  exit,  the  total  pressure  loss  across  the  bow  shock  becomes  quite  large. 
Therefore,  if  high  stagnation  pressure  is  required  on  the  model,  it  will  not  be  possible 
to  test  downstream  of  the  Mach  rhombus. 

The  stagnation  point  velocity  gradient  and  shock  standoff  distance  are  presented  in 
Figures  44  and  45.  The  greatest  sensitivity'  to  body  size  is  seen  to  be  immediately 
downstream  of  the  Mach  rhombus.  This  is  precisely  the  location  where  the  gradient 
in  flow  angle  is  most  severe.  (See  Figures  3  and  6.)  Therefore,  the  worst  possible 
simulation  of  the  flight  environment  (i.c.,  uniform  flow)  occurs  immediately  down¬ 
stream  of  the  Mach  rhombus. 
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The  results  of  the  above  calculations  indicated  that  the  Mach  No.  gradients  and  the  flow 
angle  gradient  in  the  z  direction  exerted  only  a  second  order  influence  on  the 
solution.  For  the  range  of  body  sizes  considered  »  1^,  the  flow  angle  gradient 

normal  to  the  axis  was  the  dominating  factor.  Examination  of  Figures  5  and  G  indicate 
that  the  flow  angle  is  nearly  linear  in  r  out  to  the  location  of  the  plume  barrel  shock. 
Furthermore,  in  the  vicinity  of  the  flow  field  axis,  the  influence  of  the  nonuniform 
free  stream  is  exerted  only  through  the  normal  component  of  the  flow  angle  gradient 
(i.e. ,  for  a  nonuniform  external  flow,  the  stagnation  point  Eq  (5)  and  (G)  contain  only 

d  <t>  \ 

the  additional  term,  j . 

The  above  considerations  suggest  that  the  influence  of  the  nonuniform  free  stream 

^  4> 

might  be  characterized  by  a  single  parameter,  the  value  of  - —  on  the  flow  field 

Sr 

%  4> 

axis.  For  a  body  inside  the  Mach  rhombus,  the  appropriate  value  would  be  - — 

Or 

evaluated  at  the  interface  between  uniform  and  nonuniform  flow.  In  order  to  test  this 

b4> 

hypothesis,  several  solutions  were  generated  for  constant  Mach  No. ,  -  s  0,  and  a 

dz 

d4> 

constant  value  of  .  The  approximate  solution  is  indicated  by  the  triangles  in 
or 

Figure  43.  The  error  in  the  pressure  distribution  is  approximately  10  to  15  percent. 
(For  smaller  values  of  —  ,  the  error  would  decrease.)  The  error  in  the  shock 


standoff  distance  and  stagnation  point  velocity  gradient  was  about  five  percent. 


It  appears,  therefore,  that  the  description  of  the  plume  flow  field  may  be  greatly 
simplified  with  only  a  small  sacrifice  in  the  accuracy  of  the  blunt  body  solution. 
Thus,  instead  of  specifying  the  entire  plume  flow  field,  only  the  values  along  the 
plume  axis  are  required.  Since  the  plume  axis  is  a  streamline,  the  flow  angle 
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gradient  may  be  expressed  in  terms  of  the  Mach  No,  gradient  along  the  axis.  That  is, 
for  an  axisymmetric,  inviscid,  irrotational  flow  field,  the  following  expression  is 
valid  on  the  flow  field  axis: 

d<t>  _  M2-! _  dM 

dr  2M  (l  +  —■  M2) 

In  summary',  if  it  is  assumed  that  the  effect  of  the  nonuniform  free  stream  may  be 

characterized  solely  by  the  value  of  — —  on  the  plume  axis,  then  the  description 

or 

of  the  plume  flow’  field  is  reduced  to  specifying  the  Mach  No,  and  its  gradient  on  the  axis. 


c 


20 


SECTION  G 


SUMMARY  AND  CONCLUSIONS 


It  has  been  shown  that  a  nonuniform,  expanding  flow  field  can  significantly  affect  all 
the  characteristics  of  a  blunt  body  immersed  in  that  flow  field.  In  particular,  solutions 
have  been  generated  for  a  sphere  and  circular  disk  in  flows  where  the  free  stream 
divergence  is  of  the  same  order  as  may  be  encountered  in  underexpanded  arc  jets  or 
rocket  engines. 

In  light  of  the  similarity  between  spherical  source  flows  and  underexpanded  exhaust 
plumes,  an  extensive  set  of  calculations  has  been  carried  out  over  a  range  of  body  sizes 
to  source  distances.  It  is  felt  that  by  suitably  locating  the  "effective"  origin  of  the  flow, 
these  solutions  may  be  used  to  provide  reasonable  estimates  of  the  blunt  body  character¬ 
istics  for  arbitrary  nozzle  geometries  and  operating  conditions.  Since  the  Mach 
gradient  for  the  spherical  source  flows  is  quite  weak,  reasonable  simulation  should  be 
expected  only  for  exhaust  plumes  where  the  flow  angle  gradient  is  the  dominating 
factor. 

It  has  been  demonstrated  that  there  exists  the  possibility  of  the  sonic  point  moving  off 
a  sharp  corner  onto  the  face  of  a  blunt  body,  even  for  bodies  with  zero  surface  curva¬ 
ture.  The  results  of  the  calculations  for  an  infinite  circular  disk  may  be  used  to  assess 
the  probability  of  this  occurring. 

The  description  of  the  nonuniform  free  stream  may  itself  pose  considerable  difficulties, 
particularly  if  it  is  a  rotational  inviscid  flow  or  a  viscous  How.  Even  the  irrotational 
exhaust  plume  requires  a  numerical  solution  to  the  equations  of  motion.  When  the 
definition  of  the  external  flow  is  in  terms  of  a  tabular  array  in  two  space  variables, 
considerable  effort  may  be  required  to  evaluate  all  the  derivatives  of  the  free  stream 
properties.  Furthermore,  the  blunt  body  solutions  proved  to  be  quite  sensitive  to 
small  errors  in  interpolating  for  the  free  stream  properties.  Therefore,  the  effect 
upon  the  blunt  body  solutions,  of  various  simplifications  for  the  free  stream 
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description,  has  been  examined.  Solutions  were  generated  for  a  sphere,  where  the 
Mach  No.  and  flow  angle  variation  were  correctly  accounted  for  but  the  derivatives  of 
Mach  No.  and  flow  angle  were  ignored.  The  results  indicated  that  these  derivatives  could 
be  ignored  without  introducing  large  errors.  The  simplification  of  the  free  stream 
definition  was  carried  a  step  further  for  a  few  of  the  exhaust  plume  calculations,  where 
only  the  radial  variation  of  flow  angle  was  accounted  for.  Again,  the  error  introduced 
by  this  approximation  was  relatively  small. 

In  order  to  obtain  high  stagnation  pressures  on  the  model,  most  ablation  tests  are 
conducted  at  low  supersonic  Mach  No.  Usually  the  nozzle  is  contoured  to  produce  a  region 
of  uniform  flow  downstream  of  the  exit  plane.  Since  the  Mach  rhombus  is  limited  in 
size,  it  imposes  a  size  constraint  on  the  test  model  if  uniform  flow  is  desired.  In 
order  to  examine  the  effect  of  body  size  on  the  simulated  flight  environment,  computa¬ 
tions  wei’e  performed  for  a  sphere  in  a  Mach  2  exhaust  plume.  It  has  been  shown  that 
the  maximum  body  size,  for  uniform  flow  simulation  in  the  Mach  rhombus,  is  about 
one  half  the  nozzle  exit  radius..  In  order  to  obtain  a  uniform  flow  pressure  distribution 
on  a  body  which  is  equal  in  size  to  the  nozzle  exit  radius,  it  is  necessary  to  test  four 
nozzle  radii  downstream  of  the  exit  plane.  The  region  immediately  downstream  of  the 
Mach  rhombus  has  been  identified  as  the  worst  possible  test  location  for  the  simulation 
of  uniform  flow. 

Due  to  the  difficulty  and  expense  of  obtaining  data  in  high  temperature  flows,  there 
does  not  exist  any  comprehensive  set  of  measurements  with  which  the  present  calcula¬ 
tions  can  be  compared.  Furthermore,  in  spite  of  the  large  amount  of  impingement 
testing  in  exhaust  plumes,  there  also  does  not  appear  to  be  any  cold  flow  data  which 
arc  applicable  to  the  present  problem.  This  points  up  the  need  for  additional  experi¬ 
mental  work  to  define  the  capabilities  and  limitations  of  testing  in  underexpanded 
exhaust  plumes. 
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APPENDIX  1 


INTEGRAL  RELATIONS  EOR  A  NONUNIFORM  FREE  STREAM 


It  is  the  purpose  of  this  section  to  develop  the  method  of  integral  relations  for  a  non- 
uniform  fi'ce  stream. 

The  three  basic  equations  are  listed  in  Section  3  as  Eq  1-3.  The  dependent  variables 
are  the  shock  standoff  distance,  8  ,  the  surface  velocity,  ,  and  the  shock  inclina¬ 
tion  angle,  y3.  The  independent  variable  is  distance  measured  along  the  body 
surface,  s  . 

The  oblique  shock  relations  are  required  to  evaluate  the  properties  behind  the  shock. 
These  may  be  written  as  follows: 

Pj  (7  +  1)  Sin“  (p-4>) 

(7-1)  M?,  Sin2  (£-<*>)  +  2 

2TM2  Sin2  (p- «f>)  -  (7  -  1) 
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Refer  to  Figure  1  for  sign  conventions  and  geometry. 

In  order  to  evaluate  the  derivatives  of  the  shock  properties  with  respect  to  s,  the 
partial  derivatives  of  the  properties  belli nd  the  shock  must  be  obtained. 
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The  only  remaining  derivatives  to  evaluate  are  those  involving  the  free  stream  density, 
velocity  and  dynamic  pressure: 


p  v  =  yp  iv  r 

CO  CO  CO  3 


-  y 
y-i 


_  /  (y-wi 

Vmax  ^  2  +  (T-l)M2 

CO 


V 

max 


1  + 


y-i 


-1 
-i  y-i 


M 


yPr 


(y-D  a 


Tec 


_d_ 

ds 


2m  (y-i) 

CO 


[2  +  (y-l)M^ 


dM 

a 

ds 


t 


I; 


I. 


29 


Thus,  given  the  free  stream  Mach  No.,  flow  angle,  total  pressure  and  total  enthalpy,  and 
their  gradients,  all  the  above  expressions  may  be  computed.  It  is  noted  that  if  the  total 
pressure  and  enthalpy  are  constant,  then  the  free  stream  thermodynamic  state  properties 
may  be  normalized  with  their  respective  total  values.  If  that  is  the  case,  then  only  the 
Mach  No.  and  flow  angle  are  needed  to  completely  define  the  free  stream  flow  field. 
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APPENDIX  2 


STAGNATION  POINT  BEHAVIOR  OK  THE  DIFFERENTIAL  EQUATIONS 


As  mentioned  earlier,  Eq  1  and  2  are  indeterminate  at  the  stagnation  point.  To 
evaluate  the  initial  derivatives  for  0  and  u^,  L'hopitals  rule  may  be  applied  to  the 
equations. 

At  the  stagnation  point,  the  initial  values  may  be  expressed  as  follows- 
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From  the  results  of  Appendix  1 ,  we  get: 
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Substitution  of  those  expressions  into  Eq  2  yields  the  following: 
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Again,  from  Appendix  1,  we  gel 
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Substitution  of  this  c>qn'cssion  and  Eq  2-1  into  Eq  1  yields  the  following  expression  for 
the  stagnation  point  velocity: 
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It  is  interesting  to  note  that  ^  cancels  out  upon  substitution  of  eq.  2-1  into  Eq  1. 
Therefore,  the  stagnation  point  velocity  gradient  has  the  same  form  for  cither  uniform 
or  non  uniform  external  flow. 
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LIST  OF  SYMBOLS 


s 

n 

6 

R 

e 

P 

<t> 

r 

2 

u 

V 

p 

p 

y 

(o 

M 

K 

T 


Distance1  measured  along  body  surface. 

Distance  measured  normal  to  body  surface. 

Shock  standoff  distance  -  measured  normal  to  body  surface. 

Radius  of  sphere  or  circular  disk. 

Body  surface  angle. 

Shock  angle  measured  to  the  body  axis. 

Flow  inclination  angle  between  free  stream  velocity  vector  and  flow  field  axis. 
Distance  measured  normal  to  flow  field  or  body  aisis. 

Distance  measured  parallel  to  flow  field  or  body  axis. 

Velocity  component  parallel  to  body  surface. 

Velocity  component  normal  to  body  surface. 

Static  pressure. 

Static  density. 

Ratio  of  specific  heats. 

Distance  measured  from  a  spherical  source  to  the  stagnation  point  of  a 
body  in  the  source  flow  field. 

Mach  number. 

Local  curvature  of  body  surface. 

1 

Pj  (0) 

P \ (0) 
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V  -  Maximum  thermodynamic  limiting  velocity  of  the  free  stream, 

max  h 

r  -  Nozzle  exit  radius 

1  e 

Ht  -  Total  enthalpy 

PT  -  Total  pressure 

SUBSCRIPTS: 

®  -  Free  stream  flow  conditions. 

0  -  Body  surface  conditions  or  conditions  at  the  stagnation  point  location 

of  a  body  in  a  spherical  source  flow. 

1  -  Conditions  immediately  behind  the  shock  wave. 

u  -  Uniform  flow  conditions. 


SUPERSCRIPTS: 


* 


Sonic  point  conditions  -  either  on  the  body  surface  or  on  the  bow  shock. 


Refers  to  solutions  for  which  the  free  stream  derivatives  of  mach  and 
flow  angle  were  ignored,  that  is; 


do  _ 
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Table  1 

SUMMARY  OF  SOLUTIONS  FOR  A  SPHERE  IN  A  SPHERICAL  SOURCE  FLOW 


Tabic  3 


SUMMARY  OF  SOLUTIONS  FOR  AN  INFINITE  CIRCULAR 
DISK  IN  A  SPHERICAL  SOURCE  FLOW 


■ 

STAGNATION  POINT 

SONIC  POINT  COORDINATES 

BBHBMM 

PARAMETERS 

BODY 

SHOCK  WAVE 

8 

du 

0  o 

S* 

S* 

6 * 

y 

Mo 

V 

V  ds 

max 

!  0 

(0 

1.20 

2.0 

0. 1706 

0.4745 

0.  868 

0.558 

0.220 

3.0 

0. 1141 

0. 560G 

0.  759 

0.  472 

0. 141 

5.  0 

0.0801 

0.6555 

0.  671 

0.  390 

0.  096 

10.0 

0. 0646 

0.  7224 

0.  626 

0.  348 

0.  076 

20.0 

0.0606 

0.  7436 

0.  614 

0.  336 

0.  071 

1.4 

2.0 

0. 1910 

0.  6301 

0.  865 

0.575 

0.242 

3.0 

0. 1410 

0.  7123 

0.760 

0.530 

0.176 

5.0 

0. 1119 

0.  7880 

0.688 

0.4S1 

0. 137 

10.0 

0.0991 

0. 8326 

0.654 

0.455 

0.  120 

20.0 

0.0958 

0.  8454 

0.644 

0.447 

0.116 

1. 6667 

2.0 

0.2131 

0.7584 

0.  880 

0.583 

0.  263 

3.0 

0. 1681 

0.  8294 

0.781 

0.578 

0.210 

5.0 

0. 1428 

0.8877 

0.722 

0.559 

0. 179 

10.0 

0. 1318 

0.918S 

0.693 

0.  544 

0. 164 

20.0 

0. 1291 

0.9274 

0.687 

0.  539 

0.  161 

Table  4 


SUMMARY  OF  SOLUTIONS  FOR  A  SPHERE  IN  A  MACH  2  EXHAUST  PLUME  7  1.40 


STAGNATION  POINT 
PARAMETERS 

SONIC  POINT  COORDINA'I  ES 

BODY 

SHOCK  WAVE 

2 
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-R. 

r 

e 

fi 

R 

R  du0 

^  max  ds 

e * 

JLT 

R 

e* 

1.0 

1 

0.  302 

0.  435 

41. 1 

0.593 

38.  3 

0.  353 

0.443 

44.5 

0.471 

48.5 

■ 

0.290 

0.515 

55.  0 

0.299 

02.  6 

0.236 

0.  C01 

1.50 

0.50 

0.  280 

0.520 

1.00 

0.225 

0.  029 

2.0 

0.0 

0.208 

0.500 

44.  2 

0.  382 

47.5 

0.50 

0. 181 

0.  703 

55.9 

0.251 

58.  4 

1.00 

0. 144 

0.  800 

3.0 

0 

0. 197 

0.579 

0.50 

0.  157 

0.  704 

1.0 

0.  128 

0.  847 

2.0 

0.092 

1.141 

4.0 

0 

0.  175 

0.  013 

47.  7 

0.  220 

57.3 

0.5 

0.144 

0.  724 

53.  3 

0.  1  SO 

02.  2 

1.0 

0.  123 

0.  838 

58.4 

0.  154 

65.  3 

2.0 

0.091 

1.095 

G5.2 

0.  114 
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Geometry  for  Blunt  Body  in  a  Nonuniform  Flow  Field 


Figure  2.  Flow  Field  Exhaust  Plume  for  a  Mach  2  Contoured  Nozzle 
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Figure  3.  Centerline  Properties  for  the  Exhaust  Plume 


Profiles  for  the  Contoured  Nozzle  Exhaust  Plume 
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Figure  22.  Comparison  of  Pressure  Distributions  for  a  Sphere  in  a  Spherical 
Source  Flow 
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Figure  23.  Comparison  of  Pressure  Distributions  on  a  Sphere  in  a  Spherical 
Source  Flow  with  Newtonian  Impact  Theory 
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Figure  25.  Sensitivity  of  the  Stagnation  Point  Velocity  Gradient  to  the  Free 
Stream  Gradient  Terms 
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Figure  ?.fi.  Pressure  Distributions  on  a  Circular  Disk  in  a  Spherical  Souj 
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Figure  29.  Shock  Detachment  for  a  Ci 
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Figure  31.  Mach  No.  Effect  Upon  the  Stagnation  Point  Velocity  Gradient  for  a 
Circular  Disk  in  a  Spherical  Source  Flow 
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Effect  of  Specific  Heat  Ratio  Upon  Shock  Standoff  Distance  for 
Circular  Disk  in  a  Spherical  Source  Flow 
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YHouitv  Gradient  for  an  Infinite  Circular  Disk  in  a 
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Figure  10.  Pressure  Distributions  on  a  Sphere  in  a  Mach  2  Exhaust  Plume 
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Figure  42.  Turbulent  Heat  Transfer  Distributions  for  a  Sphere  in  a 
Mach  2  Exhaust  Plume 
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